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Abstract 

Stellar seismology reveals some interior properties of thousands of solar-type 
stars but the solar seismic sound speed stays puzzling since a decade as it dis¬ 
agrees with the Standard Solar Model (SSM) prediction. One of the explana¬ 
tions of this disagreement may be found in the treatment of the transport of 
radiation from the solar core to the surface. As the same framework is used 
for other stars, it is important to check precisely the reliability of the interact¬ 
ing cross sections of photons with each species in order to ensure the energy 
transport for temperature T > 2 - 10^ K and density p > 0.2 g/cm^. In this 
paper, we propose a new technique to reach the domain of temperature and 
density found in the solar radiative interior. This technique called the Double 
Ablation Front (DAF) is based on a high conversion of the laser energy into 
X-rays thanks to moderated Z material irradiated by laser intensity between 
1.5 X 10^^ W/cm^ and 4 x 10^^ W/cm^. This high conversion creates, in 
addition to the electronic front a second ablation front in the moderated Z ma¬ 
terial. Between the two fronts there is a plateau of density and temperature that 
we exploit to heat a sample of iron or of oxide. The first simulations realized 
with the hydrodynamic code CHIC show that this technique allows to reach 
conditions equivalent to half the radiative zone of the Sun with high stability 
both in time and space. We examine the possibility to measure both iron and 
oxygen absorption spectra. 
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1. Introduction 

The SSM including the updated photospheric composition [T] in carbon, oxy¬ 
gen and nitrogen (CNO) disagrees with helioseismic radial profiles and neutrino 
detections El 0 a ig. This is particularly visible on the sound speed profile 
that is sensitive to the detailed internal solar composition through the opacity 
coefficients 0 and references therein]: Figure recalls the relative difference 
between the squared sound speed coming from seismology and from the SSM. 
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Figure 1: Relative differences between the squared sound speed coming from seismology 
(SoHO) and from the SSM (red). A seismic model has been calculated with the same equa¬ 
tions than SSM but it has been adjusted to respect as properly as possible the seismic results. 
Its relative differences with observations are drawn with a black continuous line, the seismic 
error bars are shown on that model. Adapted from Ref. [ 3 ] (Color online). 


This discrepancy between the observed sound speed and the sound speed 
predicted by the SSM varies largely greater than the vertical error bar [3l[5]. To 
explain this difference, three hypotheses, which could all exist simultaneously, 
have been advanced IS1I7]. One can first question the transfer of energy from 
the center of the Sun to the surface through the opacity coefficients: the atomic 
calculations in plasma conditions could underestimate the Rosseland mean opac¬ 
ities which directly drive the radiative transfer. The second idea is to put some 
doubt on the gravitational settling which could be underestimated for CNO and 
heavy elements due to an incorrect treatment of the radiative acceleration of 
elements towards the surface with, as a consequence, incorrect central abun¬ 
dances. The last hypothesis concerns the energetic balance equation: the Sun 
could produce slightly more energy (< 5 %) than it liberates at the surface and 
this additional energy would be transformed into macroscopic motions existing 
in the radiative zone [8], so the energetic balance of the SSM is incomplete as 
it does not contain any internal dynamical effects. 

Determining the origin of this discrepancy (in the opacity ingredients or 
in some limitation of the solar model) would be an important step toward a 
better understanding of the solar interior. As the two first hypotheses put 
in question the way the radiative transfer is used and as the third hypothesis 
could be difficult to verify except by a very precise neutrino detection, it is 
important to test the opacity calculations in details. So, this situation requires 
an experimental validation of the radiative transfer calculations at the conditions 
of the solar radiative zone (see Table for the solar internal conditions to be 
explored). 
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Figure 2: Relative contribution of the most important heavy element to the total Rosseland 
opacity (including H and He) for the internal conditions of the Sun, the composition of Ref. 
using OPAL opacities [9]. From Ref. [8]. 


2. An experimental challenge 

The Sun is principally constituted of hydrogen and helium which are, in al¬ 
most the whole Sun completely ionized. However, even though heavy elements 
are present only at a few percents in mass (iron represents only about 10“^ of 
the hydrogen contribution in mass fraction), they contribute significantly to the 
global opacity [8]. 

Figure represents the respective role of the main elements contributing to 
the global opacity along the solar radial profile. The elements heavier than 
^He represent only 1.4 % in mass fraction but their contribution to the opac¬ 
ity coefficient is up to 70 % at the basis of the convective zone. The largest 
contribution of iron and oxygen in the radiative zone is due to their bound- 
bound contributions. Iron is never completely ionized in the Sun: even in the 
center, it still contributes to opacity with bound-bound and bound-free transi¬ 
tions. Its contribution to the total opacity is around 20 - 30 % in most of the 
radiative zone. Oxygen, the third element in abundance, changes from fully to 
partly ionized above 0.4 Rq and plays a major role to trigger the convection 
instability around 0.7 R©. Unfortunately, these plasma properties have never 
been verified in laboratory nor the absorption energy spectra. Several problems 
have made such measurements difficult. To properly evaluate the opacity of the 
considered elements, one has to reproduce their charge state distribution and 
the free-electron density Ng at the targeted conditions. As shown in Table 
the free-electron densities in the radiative zone are 10^^ — 10^"^ cm“^. Then, 
one needs to determine monochromatic single-element opacity in well diagnosed 
plasmas with good uniformity of temperature and density. Figure shows the 
charge state distributions of iron and oxygen at the conditions of Table I. As 
iron is partially ionized, with a great number of bound electrons, the calculation 
of its opacity is challenging. 
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Table 1: Summary of solar conditions found in Ref. m and in MESA solar model in 


Solar radius (r/RQ) 

T(eV) 

p (g.cin-3) 

Ne (cm“3) 

0.5 

340 

1.36 

8 X 10^3 

0.6 

270 

0.50 

2.5 X 10^3 

0.7 

200 

0.21 

1 X 10^3 




Figure 3: Charge state distributions of iron (left) and oxygen (right) at 0.5 (...), 0.6 (—) and 
0.7 (full line) Rq obtained with FLYCHK | 12| . 


In the solar radiative zone, the radiation transport is described by a diffusion 
approximation using the Rosseland mean opacity kr: 


J_ ^ 

kr 

where B is the Planck function, T is the temperature and the spectral 

opacity. The shape of the weighting function dB/dT is represented on Figure 
1^ for the same three conditions. This Rosseland mean ponderation directly 
determines the spectral range of interest: the maximum is around to hz//kT :^4, 
it determines the part of the spectrum which contributes the most to opacity. 
So, for the solar conditions of Table[^ the range of interest is principally between 
500 and 2500 eV as previously shown by m- This gives precious information 
for the qualification of diagnostics. 


3. State of the art 

A first experiment was performed at the Z-pinch facility of Sandia National 
Laboratory at Te = 156 ± 6 eV and Ne = 6.9 ± 1.7 x 10^^ cm“^ where the 
transmission of a mixed Mg and Fe plasma was measured [14]. Recently, the 
same group performed another experiment with the same principle also on a 
mixed Mg and Fe sample and they reached Tg = 196 ± 6 eV and Ng = 3.8 ± 
0.8 X 10^^ cm' -3 [IS]. These conditions are not so far from the solar ones but 
there is a clear discrepancy between experiment and all the opacity calculations 
up to now. 
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Figure 4: Rosseland mean ponderation at the three different temperatures. 


Other experiments involving lasers and compared to theoretical opacities 
[HiniiiH] were limited to Tg < 100 eV and density p < 10 ^ g/cm which 
is far from solar conditions. To heat material at Te>200 eV at relatively high 
density, there are different methods with laser beams: with a shock (giving Tg of 
few tens eV and p ^ 3 - 4 times the solid density, with small gradients) [19] , with 
a thermal wave (Tg ^ few tens eV and density around the solid density, with 
high gradients) [19] and finally with short pulse lasers (Tg ^ 600 - 800 eV and 
p > solid density) [20l|2T]. There is also another technique, proposed in Ref. 
[22] . using X-rays conversion in holhraum, which is presently in development. 
They hope to obtain Tg of some hundreds eV and p around 0.1 the solid density. 

We propose in this paper another approach, called the Double Ablation Front 
(DAF) to limit the gradients in the foil. This method is promising for measuring 
iron absorption spectrum in solar conditions. We explain in the next sections 
the interest of this technique and show some characteristics of the design that 
we apply here to an oxide measurement. This DAF has been experimentally 
put in evidence by [23] and extensively studied for fusion by [24] and [25]. 

4. CHIC code 

We have used the CHIC code (Code d’Hydrodynamique et d’Implosion du 
CELIA) [26] for our simulations. This hydrodynamic code is Lagrangian and 
dedicated to ICF calculations. It operates in ID and 2D axially symmetric 
geometries and includes three-dimensional ray tracing for laser beam propaga¬ 
tion. The laser energy release is modeled by inverse bremsstrahlung absorption. 
The electronic thermal transport is described in the classical Spitzer-Harm ap¬ 
proximation, with a flux limitation of 0.06 [23 [28]. In our calculations, the 
wavelength of the laser beams is set to 0.351 pm corresponding to both LMJ 
and NIF beams, the diameter of the spot being 400 pm. The pulse duration 
is 2.2 ns with a linear ramp of 0.2 ns at the beginning and at the end of the 
pulse. The SESAME table [29] is used for the equation of state and we con¬ 
verge in number of groups for the treatment of the radiative transport through 
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the diffusion approximation. In the following simulations, we did not take into 
account any magnetic field generation nor a non local conduction. 


5. Double Ablation Front 

The DAF approach is based on an enhancement of the radiative effects in 
the target, due to a high conversion between the laser energy and X-rays in 
a moderated Z material. The laser when deposing its energy near the critical 



Figure 5: Schematic profiles of electronic temperature Te, density p and opacity hiR at a given 
time in a layer of moderated Z ablator. One sees different regions: the corona (Ne< Nc] Te 
around 3 keV and p ~ 0.02 g/cc); the conversion zone (Te around 1 keV and low density); 
the reemission zone (Te of several hundreds eV and p ~ 1 g/cc), heated by the X-rays | 30| : 
the shock wave and the non-perturbed solid. Two ablation fronts appear, one due to electrons 
(electronic front), the other due to photons (radiative front) (Color online). 


density Nc strongly heats this region. X-rays are emitted, through several pro¬ 
cesses (principally Bremsstrahlung emission (free-free transitions), emission by 
electron-ion recombination (bound-free transitions) and line emission (bound- 
bound transitions), dominant in the case of high Z materials). These X-rays 
are isotropically emitted: the radiation propagating toward the target in high- 
density region (^ solid density) and relatively low temperature (^ hundreds 
eV) will be highly absorbed because of the high opacity of this region, creating 
a radiative ablation front [24l|25]. The higher the atomic number Z is, the more 
important will be the X-rays emission (Bremsstrahlung emission proportional 
to 7?). At the same time, the electrons also transport the energy through the 
target, beyond the critical density. The electronic temperature decreasing in 
the denser region, they deposit their energy, creating a thermal ablation front 
(see Figure [^. Between the two fronts, there is a ’’plateau” of density which 
can be exploited as presented in the following section. Because of the difference 
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Figure 6: Left : One-side irradiation target, composed of one layer of ablator (green), one 
layer of the sample of interest (yellow) and one layer of tamper (red). Right : symmetrical 
irradiation target, composed of two layers of ablators (green) tampering the sample of interest 
(yellow) (Color online). 


of velocities of the two fronts, this plateau region extends with time, ensuring 
therefore its stability for diagnosis during an experiment. The diffusion approx¬ 
imation is satisfying in this ablation region because of its high opacity which 
reduces drastically the photon mean free path. This hypothesis is no longer valid 
in the low optical depth region of the corona but it has already been shown that 
the diffusion approximation yields to quantitatively good results [31]. 

This multi-ablation structure was studied experimentally by [23] in 2004. They 
observed the evolution of a DAT structure inside a CHBr target and showed 
that this structure leads to the suppression of the Rayleigh-Taylor instability, 
because of the larger ablation velocity and the longer density scale length. This 
structure was also produced more recently at the OMEGA laser by [32] . 

6. Experimental aspects 

We performed simulations on three-layer planar targets with two different 
designs: the first type is constituted of one layer of ablator, creating the DAF 
structure, followed by the sample of interest, tampered by a layer of CH and 
irradiated by only one laser beam; the second design is constituted of two layers 
of ablator, tampering the sample of interest and attacked symmetrically by two 
laser beams (Figure]^. 

The objective is to measure the transmission of X-rays through the sample 
that is related to the opacity by: 

where p is the density, n the spectral opacity and r the thickness of the probed 
plasma. 

We extract from the simulation the mean density Pmean^ temperature 
and free electron density A^emean ^ defined by: 

I Pdm jTedm J Nedm 

/ dm / dm / dm 

with their gradients on the plateau and in the iron foil, dm is the mass element. 
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7. Measurements of iron absorption spectra 

Choice of the ablator. We have investigated different materials to find the best 
one for heating iron and oxygen at the solar conditions. We performed sim¬ 
ulations with the same target design, modifying only the type of the ablator. 
We use a laser intensity liaser=^ x 10^^ W/cm^, irradiating three-layer targets: 
ablator (10 jam) / Fe (1 jam) / CH (10 jam) and one follows the evolution of 
the mean density, the mean temperature and the mean free electron density 
in the sample (see Table . We note that high Z ablators (Z > 20) lead to 


Table 2: Temperature, density and free electrons density in an iron sample during the plateau 
region obtained with different ablators and a laser irradiation of 1 x 10^^ W/cm^. 


Ablator 

Z 

T(eV) 

P (g/cm-^) 

Ne (cm ^) 

Titanium 

22 

150 

0.50 

0.6 X 10^^ 

Copper 

29 

130 

0.25 

1 X 10^^ 

Silicon 

14 

155 

0.90 

1.4 X 10^3 

Sapphire 

-14.8 

140 

1.16 

1.7 X 1023 

Quartz 

10 

153 

1.14 

1.6 X 1023 

Aluminium 

13 

155 

0.85 

1.3 X 1023 


important radiative effects (high conversion between the laser energy and the 
X-rays). Consequently, the sample is strongly pre-heated and strongly expands. 
When the shock goes through the sample, the increase of the density is hence 
less important in comparison with a moderated Z ablator. The density reached 
during the plateau is then smaller (in the case of titanium (Z = 22): 0.5 g/cm^; 
in the case of quartz (Z = 10): 1.14 g/cm^). 

To discriminate between the different ablators, the spectra of the transmitted 
photons have to be considered. The objective is to avoid perturbation in the 
energy range of interest. Figure shows the comparison of the opacity of the 
different materials during the plateau region in the sample. Silicon is the best 
adapted ablator as it does not perturbed the region between 500 and 1500 eV 
where the important structures that we would like to study are present. So the 
following simulations use silicon as an ablator. 

Irradiation with one beam. For a Si (8 /am) / Fe (0.1 /am) / CH (7 fim) target, 
irradiated by a laser intensity of 1.5 x 10^^ W/cm^, we obtain 0.75 < Pmean < 
1 g/cm^, 160 < Tmean < 180 and 1.1 < < 1.5 X 10^^ cm“^ in the 

sample on a plateau of 0.7 ns. The spatial gradients are no more than 8 % in 
the sample. 

Irradiation with two beams. In the case of a symmetrical irradiation by two 
beams of same intensity 1.5 x 10^^ W/cm^ arriving on a target of Si (7 jam) 
/ Fe (0.1 jam) / Si (7 /im), we got a density of 1.2 < Pmean <1-5 g/cm^, 198 
< Te^g^^ < 228 eV and between 2.2 < < 2.50 x 10^^ cm“^ with gradients 










Figure 7: Opacity of the considered ablators, in typical conditions of density and temperature 
of the simulated experiment (at the transmission measurement time) (Color online). 



Figure 8: Mean density, temperature and free electron density profiles obtained with a sym¬ 
metrical laser irradiation of intensity 1.5 X 10^^ W/cm^ on Si (7 /rm) / Fe (0.1 /rm) / Si (7 
p,m) (Color online). 


lower than 5 % on 0.7 ns in the iron sample (see Figure [^. 
Figure shows the transmission of the whole target. The choice of silicon 
as an ablator is validated as the structure of interest of the iron spectrum is 
clearly separated from silicon ones at these temperatures, which would allow a 
clear measurement. The charge state distribution at these conditions is pretty 
close from the solar one (around 20 % differences on the fraction of the relevant 
ions (Fe XVI to Fe XIX) with the charge state distribution at 0.7 R©), providing 
therefore a meaningful test of opacity calculations at these conditions. 

The two previous simulations were performed with a laser intensity close 
to the first LMJ configuration. In the future or with NIF, one can expect to 
reach higher intensities. So a third simulation, a double irradiation with a laser 
intensity of 4 x 10^^ W/cm^ with the same target as the previous case, has 
been also computed. The results of the simulation give between 2.0 and 2.3 
g/cm^ in mean density, between 265 and 290 eV in temperature and between 
3.7 and 4.1 x 10^^ cm“^ for All spatial gradients were under 3 %, due 

to the huge compression induced by the lasers on each side of the target. The 
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Figure 9: Transmission of iron and silicon at the measurement time obtained with a symmet¬ 
rical laser irradiation of intensity 1.5 xlO^^ W/cm^ on the Si (7 fim) / Fe (0.1 fim) / Si (7 
/rm) target (Color online). 


charge state distribution exhibits around 4 % of difference on the main ion (Fe 
XVIII) and around 15 - 20 % for the others with the charge state distribution of 
iron at 0.6 R©. So, with even higher intensities, one can expect to reach higher 
temperatures and densities, with very high stability in space and time in the 
sample and then might reproduce conditions corresponding to the middle of the 
radiative zone. 


8. Measurements of oxgen absorption spectra 

Oxygen is a particularly interesting element to study for solar application. 
Measuring the opacity of pure oxygen is really complicated at these conditions 
of high density and temperature : one could use a gas but then the reached 
density would not be high enough. That’s why we propose to use oxides. We 
have performed simulations with the same design as the one used for iron, with 
hematite as a sample (Fe 203 ): Si (7 fim) / Fe 203 (0.1 fim) / Si (7 /rm). We 
obtained between 0.87 and 1.08 g/cm^ in mean density, between 211 and 242 eV 
for the mean temperature and between 1.95 and 2.3 x 10^^ cm“^ for the free 
electron density (see Figure pR]). The transmission of oxygen, iron and silicon 
at these conditions is represented on Figure [TT| We are presently studying the 
relation between the transmission of oxygen and iron and the transmission of 
Fe203. 

We are also performing simulations on other types of oxides to optimize the 
oxygen opacity measurement. 
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Figure 10: Mean density, temperature and free electron density profiles obtained with a 
symmetrical laser irradiation of intensity 1.5 X 10^^ W/cm^ on Si (7 p,m) / Fe 203 (0.1 p,m) 
/ Si (7 /rm) (Color online). 
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Figure 11: Transmission of oxygen, iron and silicon in the conditions of the plateau obtained 
with a symmetrical laser irradiation of intensity 1.5 xlO^^ W/cm^ on the Si (7 p,m) / Fe 2 03 
(0.1 fim) / Si (7 fim) target (Color online). 


9. Conclusion 

Discrepancies between the seismic observations and the prediction of the 
solar standard model requires an experimental validation of opacity calculations. 
The DAF approach creates conditions equivalent in charge state distribution and 
free electron density to the conditions of the radiative zone (Te around 200 - 
300 eV and Ne ^ few 10^^ cm“^), with high stability of the density and the 
temperature both in time and space (gradients smaller than 10 %), convenient 
for LTE measurements and to check the opacity calculations for specific elements 
or mixture. With 1.5 x 10^^ W/cm^ (corresponding to first LMJ-PETAL facility 
configuration), it is possible to reach interesting astrophysical conditions with 
high stability. With the addition of quads on LMJ-PETAL or with NIE, one 
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will increase temperature and density to begin to look at plasma effects. 

After the first experimental validation of this concept, the check of the line 
broadening will require high-resolution spectrometer for a definitive validation 
of the detailed calculations. This work will benefit to inertial fusion by limiting 
instability development in different experimental studies. 
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